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ABSTRACT 

A  model  for  bistatic  acoustic  scattering  at  high  frequencies  (10-100  kHz)  from 
elastic  seafloors  is  developed  by  combining  the  elastic  small-slope  roughness  scattering 
model  of  Yang  and  Broschat  with  the  elastic  volume  scattering  model  of  Ivakin. 
The  combined  model  is  applicable  to  a  variety  of  seafloors,  including  clays,  silts, 
sands,  gravel,  and  rock.  Although  elastic  effects  are  negligible  in  sands  and  finer 
sediments,  they  are  essential  in  treating  reflection  and  scattering  by  rock,  which 
cannot  be  usefully  represented  by  a  fluid  model.  Use  of  the  small-slope  approximation 
makes  it  possible  to  avoid  the  cumbersome  interpolation  between  the  Kirchhoff  and 
perturbation  approximations  employed  in  the  previous  APL-UW  scattering  models 
and  allows  larger  levels  of  roughness  than  are  permitted  with  conventional  scattering 
approximations.  The  model  results  agree  with  those  of  the  APL-UW  backscattering 
model  and  with  those  of  the  APL-UW  bistatic  model  for  clays,  silts,  and  sands.  The 
inclusion  of  elastic  effects  attributes  the  high  levels  of  scattering  observed  for  cobble 
and  rock  seafloors  to  a  combination  of  roughness  and  volume  scattering,  and  replaces 
the  empirical  treatment  of  cobble  and  rock  seafloors  previously  used  in  the  APL-UW 
backscattering  model  by  a  physical  one.  For  these  cases,  the  new  model  sometimes 
produces  unphysical  results  for  scattering  near  the  specular  direction,  indicating  the 
need  for  further  work  on  this  problem. 
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1.  INTRODUCTION 

The  APL-UW  high-frequency  acoustic  models  document  [1]  contains  models  for 
seafloor  reflection,  backscattering,  and  bistatic  scattering.  These  models  are  eclectic 
combinations  of  models  found  in  the  literature,  and  all  treat  the  seafloor  material 
as  a  fluid.  The  reflection  model  uses  the  flat-interface  Rayleigh  reflection  coeflacient, 
while  backscattering  is  described  by  a  combination  of  four  methods:  the  Kirchhoff 
approximation,  the  composite  roughness  approximation,  an  empirical  formula  for 
large  roughness,  and  a  semiempirical  treatment  of  sediment  volume  scattering.  The 
bistatic  model  [1]  uses  a  combination  of  three  models:  perturbation  theory  and  the 
Kirchhoff  approximation  for  roughness  scattering,  and  volume  perturbation  theory 
for  sediment  volume  scattering.  The  backscattering  and  bistatic  models  have  been 
successful  in  describing  experimental  data  from  a  wide  class  of  seafloor  types  [2- 
8],  but  certain  shortcomings  are  evident.  The  most  obvious  problem  is  that  the 
fluid  approximation  is  very  poor  for  rocky  seafloors,  where  shear-wave  effects  become 
important  [9-11].  A  related  difficulty  is  that  the  Rayleigh  reflection  coefficient  for  the 
two-fluid  interface  is  known  to  be  a  bad  approximation  for  rock  [12].  As  backscatter 
is  a  special  case  of  bistatic  scattering,  it  is  possible  to  compare  the  predictions  of 
the  APL-UW  backscatter  and  bistatic  models  directly.  Figure  1  shows  that  the  two 
models  yield  substantially  different  results  for  rougher  seafloors,  for  which  the  APL- 
UW  backscattering  model  uses  an  empirical  fit  and  the  bistatic  model  does  not.  Figure 
2  shows  that  the  models  are  very  similar  for  sandy  and  silty  seafloors.  The  differences 
at  small  grazing  angles  are  due  to  the  use  of  composite-roughness  averaging  in  the 
backscatter  model. 

This  report  describes  a  single  model  that  treats  both  backscattering  and  bistatic 
scattering  for  a  wide  range  of  seafloor  types.  This  model  combines  two  recent  devel¬ 
opments  in  the  theory  of  scattering  by  random  elastic  media.  These  newer  methods 
promise  to  alleviate  some,  if  not  all,  of  the  problems  just  mentioned  while  providing 
a  more  streamlined  formalism  that  is  less  ad  hoc  and  empirical.  This  model  is  essen¬ 
tially  equivalent  to  the  older  models  in  the  regime  where  the  latter  are  believed  to  be 
physically  meaningful.  Success  outside  this  regime  cannot  be  judged  at  present  owing 
to  a  lack  of  suitable  acoustic  and  physical  data;  however,  there  is  no  doubt  that  the 
physical  basis  of  the  new  model  is  stronger  than  that  of  the  old. 

For  scattering  by  seafloor  roughness,  the  new  model  uses  the  small-slope  formal¬ 
ism  of  Voronovich  [13]  as  adapted  to  elastic  seafloors  by  Yang  and  Broschat  [9]  and 
Wurmser  [14] .  For  seafloor  volume  scattering,  the  elastic  perturbation  approximation 
developed  by  Ivakin  [16]  is  used  in  the  form  outlined  and  extended  by  Ivakin  and 
Jackson  [11, 12].  These  approximations  will  be  applied  to  a  wide  variety  of  seafloor 
types  and  to  all  angles,  avoiding  the  interpolation  required  in  the  older  models  men¬ 
tioned  previously.  It  should  be  realized  that  the  regimes  of  accuracy  of  these  newer 
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Comparison  of  TR  9407  Scattering  Models 
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Figure  1:  Comparison  of  backscattering  strengths  predicted  by  the  APL-UW 
backscattering  and  bistatic  scattering  models  for  the  extreme  seafloor  types  cate¬ 
gorized  as  rough  rock,  rock,  cobble,  and  sandy  gravel. 


Comparison  of  TR  9407  Scattering  Models 


Grazing  Angle  (deg) 


Figure  2:  Comparison  of  backscattering  strengths  predicted  by  the  APL-UW 
backscattering  and  bistatic  scattering  models  for  the  common  seafloor  types  cate¬ 
gorized  as  coarse  sand,  medium  sand,  fine  sand,  and  silt. 
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methods  have  not  yet  been  determined.  Both  numerical  calculations  using  “exact” 
codes  and  experimental  data  are  needed  to  settle  such  issues. 

The  organization  of  this  report  is  as  follows:  Section  2  defines  the  environmen¬ 
tal  parameters  that  serve  as  model  inputs.  These  adhere  as  closely  as  possible  to 
the  inputs  of  the  earlier  models,  although  the  inclusion  of  shear  effects  necessarily 
increases  the  number  of  parameters.  Section  3  gives  essential  acoustic  definitions, 
Section  4  gives  formulas  for  the  roughness  scattering  part  of  the  model,  and  Section  5 
gives  formulas  for  the  volume  scattering  part.  Numerical  illustrations  are  provided 
in  Section  6,  including  comparisons  with  older  models.  Finally,  Section  7  contains 
conclusions. 
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2.  MODEL  INPUT  PARAMETERS 

The  model  input  parameters  can  be  broken  into  three  subsets:  those  describing  the 
average  properties  of  the  seafloor  material,  those  describing  the  roughness  of  the  in¬ 
terface  between  the  water  and  the  seafloor,  and  those  describing  the  inhomogeneity 
of  the  seafloor  material.  To  provide  the  maximum  level  of  generality  allowed  by  the 
present  understanding  of  acoustic  scattering  by  the  seafloor,  the  number  of  model 
parameters  has  been  deliberately  allowed  to  be  rather  large,  with  a  total  of  16  pa¬ 
rameters.  In  most  applications,  six  of  these  parameters  can  be  assigned  simple  default 
values.  In  sands,  silts,  and  clays,  two  additional  parameters  specifying  elastic  prop¬ 
erties  can  be  set  to  fixed  defaults,  leaving  eight  parameters  for  most  applications. 
Even  this  number  is  too  large  except  where  considerable  effort  has  been  expended  to 
measure  physical  properties  of  the  seafloor.  For  more  typical  applications,  a  method 
for  obtaining  “best-guess”  parameter  values  is  outlined  at  the  end  of  this  section.  It 
uses  tables  and  equations  given  in  the  APL-UW  models  document  [1]. 

The  model  is  frequency  dependent;  thus  the  acoustic  frequency,  /,  is  also  an 
input.  In  addition,  the  model  requires  a  value  for  the  speed  of  sound,  Cyj,  in  the  water 
immediately  above  the  sediment,  although  the  results  depend  only  weakly  on  this 
parameter. 

2.1  Parameters  Describing  Average  Seafloor  Properties 

There  are  five  model  parameters  describing  the  average  properties  of  the  seafloor 
material.  It  is  assumed  that  there  are  no  gradients  in  these  average  properties,  which 
are  parameterized  by  the  following  dimensionless  ratios: 

Op  =  (sea  floor)  /  (water)  density  ratio 

Up  =  (sea  floor  p-wave)/ (water)  sound  speed  ratio 

Ut  —  (sea  floor  s-wave)  /  (water)  sound  speed  ratio 

Sp  =  (imaginary)  /  (real)  ratio  for  p- wavenumber 

St  =  (imaginary)/(real)  ratio  for  s- wavenumber. 

The  subscripts  “p”  and  “f  ’  are  used  to  denote  compressional  and  shear  waves,  respec¬ 
tively.  In  geophysics,  compressional  waves  are  designated  “p”  waves  and  shear  waves 
are  designated  “s”  waves.  The  subscript  “t”  (for  “transverse” )  is  used  here  to  denote 
shear,  as  the  subscript  “s”  is  used  to  designate  the  scattered  field.  The  parameters  Sp 
and  St  will  be  referred  to  as  “loss”  parameters  and  are  related  to  the  corresponding 
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attenuation  coefficients  ap  and  ccs  (with  units  of  dB/m)  as  follows: 


Sp 


apUpCw  In  10 
407r/ 


(1) 


.  _  atUtC^  In  10 
*  407r/ 


(2) 


It  is  essential  that  consistent  units  be  employed.  That  is,  if  ap  is  in  dB/m  and  Cyj  is  in 
m/s,  /  must  be  in  Hz.  The  parameters  Up,  Up,  and  5p  are  identical  to  the  parameters 
u,  p,  and  5  of  the  older  models  [1].  It  is  reasonable  to  set  Ut  and  5t  to  zero  (or,  to 
avoid  numerical  difficulties,  set  i/j  to  a  value  much  smaller  than  unity)  in  fine-grained 
sediments  such  as  sands,  silts,  and  clays,  where  shear  effects  are  negligible  [11]. 


2.2  Parameters  Describing  Roughness  of  Water-Seafloor  In¬ 
terface 

The  roughness  parameters  are  identical  to  those  employed  in  earlier  fluid- 
sediment  versions  of  this  model  [1,3,4]  in  which  the  two-dimensional  spectrum  of 
roughness  is  taken  to  be  of  the  form 

W(K)  =  H  .  (3) 

The  parameter  W2  is  the  “roughness  spectral  strength”  with  units  of  m^"'''^,  and  the 
parameter  72  is  the  “roughness  spectral  exponent,”  with  values  ranging  from  2  to 
4.  In  earlier  models  [1]  this  parameter  was  denoted  7.  It  is  assumed  that  seafloor 
roughness  statistics  are  isotropic  and  Gaussian.  Isotropy  is  implicit  in  the  assumed 
dependence  of  the  spectrum  on  the  magnitude,  K,  of  the  two-dimensional  wavevector, 
K.  The  spectrum  is  normalized  such  that  its  integral  over  a  given  region  of  K-space 
is  equal  to  the  roughness  variance  for  those  spectral  components.  The  spectrum  is 
double-sided;  that  is,  the  wave  vector  argument  spans  both  positive  and  negative 
values. 


2.3  Parameters  Describing  Inhomogeneity  of  Seafloor 

The  volume  inhomogeneity  parameters  are  similar  to  those  used  by  Ivakin  and 
Jackson  [11, 12].  The  density  fluctuation  spectrum  is  assumed  to  be  of  the  form 

This  is  the  spectrum  for  the  dimensionless  ratio  of  density  fluctuations  divided  by 
the  mean  density.  As  such  it  has  units  of  m^,  and  the  “volume  spectral  strength,” 
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W3,  has  units  of  The  parameter  73  is  the  “volume  spectral  exponent,”  with 

values  ranging  from  3  to  5.  Unlike  the  roughness  spectrum,  the  volume  spectrum 
is  allowed  to  be  anisotropic,  with  anisotropy  characterized  by  the  “aspect  ratio,”  a. 
This  parameter  is  the  ratio  of  the  horizontal  length  scale  of  inhomogeneity  to  the 
corresponding  vertical  scale  and  is  generally  equal  to  or  greater  than  1  [17, 18].  While 
the  volume  inhomogeneities  may  be  anisotropic  with  respect  to  directions  in  a  vertical 
plane,  they  are  assumed  to  be  isotropic  with  respect  to  directions  in  the  horizontal 
plane;  that  is,  they  are  “transversely  isotropic.”  As  a  result,  the  bistatic  interface 
scattering  cross  section  will  be  isotropic.  The  parameter  go  sets  the  largest  scale 
of  the  density  fluctuations  [11, 18].  In  the  examples  presented  here,  this  parameter 
will  be  assigned  a  frequency-dependent  default  value,  defined  in  Section  2.4.  Gaussian 
statistics  need  not  be  assumed  for  volume  fluctuations,  as  the  first-order  perturbation 
approximation  used  here  is  applicable  for  arbitrary  statistics. 

In  addition  to  fluctuations  in  sediment  density,  the  model  also  allows  fluctuations 
in  the  compressional  and  shear  wave  speeds.  These  fluctuations  are  also  normalized 
by  division  by  the  respective  means  and  are  assumed  to  have  the  same  power-law 
behavior  as  Eq.  4.  These  fluctuations  may  be  correlated  with  the  density  fluctuations, 
with  correlations  described  by  appropriate  cross  spectra.  These  assumptions  lead  to 
the  set  of  model  parameters  defined  by  the  following  equation: 

=  0(3/3' kUpp(k)  •  (^) 

In  Eq.  5,  the  dimensionless  parameters  apj}/  define  the  proportionality  between  the 
density  fluctuation  spectrum  and  the  other  volume  fluctuation  spectra, 

Wpp{k)  =  p-wave  fluctuation  spectrum 

H'tt(k)  =  s-wave  fluctuation  spectrum 

Wpp{k)  =  (p-wave)-(density)  fluctuation  cross  spectrum 

Wtp(k)  =  (s-wave)- (density)  fluctuation  cross  spectrum 

Wpt(k)  =  (p-wave)-(s-wave)  fluctuation  cross  spectrum. 

This  list  exhausts  the  possible  spectra  in  the  present  problem,  as  the  cross  spectra 
are  symmetric  in  their  subscripts;  that  is,  W(!3/3'(k)  =  The  price  of  including 

this  level  of  generality  in  the  model  is  an  additional  five  parameters,  ^pp^  ^tp) 

and  apt  (note  Upp  =  1  by  definition).  The  parameters  app  and  Ott  must  be  positive, 
while  the  parameters  Cpp,  ajp,  and  apt  are  positive  or  negative,  depending  on  whether 
the  corresponding  fluctuations  are  correlated  or  anticorrelated.  In  most  of  the  il¬ 
lustrations  presented  in  this  report,  all  five  of  these  parameters  are  set  to  zero,  as 
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in  unconsolidated  sediments  normalized  wave  speed  fluctuations  are  generally  much 
smaller  than  normalized  density  fluctuations.  This  is  a  reasonable  default  when  no 
information  is  available  on  wave  speed  fluctuations.  For  rough  rock  and  rock,  these 
parameters  are  assigned  values  adapted  from  Ref.  [11]. 

2.4  Default  Parameters  for  Sands,  Silts,  and  Clays 

This  section  will  conclude  with  suggestions  for  the  choice  of  model  parameters 
when  little  information  is  available  as  to  seafloor  type.  These  defaults  apply  to  the 
most  common  sediment  types  (sands,  silts,  and  clays).  For  extreme  seafloor  types 
(rock,  cobble,  and  gravel),  it  is  recommended  that  the  parameter  values  used  for 
illustrative  purposes  later  in  this  report  be  employed  unless  better  information  is 
available. 

As  a  first  step,  the  shear  parameters  Ut  and  6t  are  assigned  values  that  render 
shear  effects  negligible  without  introducing  singularities  into  calculations.  Values 
ut  =  0.002  and  5*  =  1  are  used  in  this  report,  although  the  only  requirement  is  that 
Ut  «  1  and  that  St  not  be  large  compared  to  unity.  The  volume  spectral  parameters 
Qpp)  (hu  O’ppi  O'tpi  and  apt  should  be  set  to  zero,  while  the  aspect  ratio,  a,  should  be 
assigned  a  value  of  unity. 

The  volume  spectrum  scale  parameter,  go,  is  inversely  proportional  to  the  length 
scale  of  the  largest  fluctuations,  and  should  be  assigned  the  small,  frequency-dependent 
value 

go  =  O.OOl/jiu  ,  (6) 

where  ky,  is  the  acoustic  wavenumber  in  water,  defined  later  in  in  Eq.  15.  This 
default  is  unphysical,  but  avoids  difficulties  that  arise  in  perturbation  theory  when 
the  volume  spectra  approach  large  values  as  the  wavevector  argument  approaches 
zero. 

For  the  six  parameters  Op,  Up,  Sp,  W2, 72  and  W3,  the  models  document  [1]  should 
be  consulted  for  default  values.  It  should  be  remembered  that  in  Ref.  [1]  Op  is  repre¬ 
sented  by  the  symbol  p  and  that  Vp  and  Sp  appear  without  subscripts.  The  volume 
inhomogeneity  exponent  is  set  to  the  value 

73  =  3  ,  (7) 

which  yields  a  volume  scattering  contribution  that  is  frequency  independent,  provided 
all  other  model  parameters  are  frequency  independent.  A  formal  difficulty  arises  with 
this  exponent  value,  as  it  yields  an  infinite  value  for  the  integral  of  the  spectrum, 
corresponding  to  an  infinite  density  ratio  variance.  Thus,  even  if  the  actual  spectrum 
follows  the  assumed  power  law  over  the  wavenumber  range  of  interest,  it  must  fall 
off  more  rapidly  than  the  third  power  of  wavenumber  at  large  wavenumbers  (short 
wavelengths).  If  we  take  this  cutoff  to  correspond  to  twice  the  acoustic  wavenumber 
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of  the  highest  model  frequency  (100  kHz),  and  take  the  lower  cutoff  to  be  the  scale 
parameter,  Eq.  6,  evaluated  at  the  lowest  model  frequency  (10  kHz),  the  approximate 
density  ratio  variance  is 

(7^  ^  u;3ln(20000)  ^  10Li;3  .  (8) 

As  the  mean  density  ratio  lies  in  the  range  1-2.5,  fluctuations  in  this  parameter  must 
be  much  less  than  unity,  leading  to  the  condition 


W3  «  0.1  . 


(9) 


This  condition  is  satisfied  by  all  the  spectral  strengths  used  in  this  report  (see  Table  1). 
The  values  for  silt  and  sand  were  obtained  using  the  default 


_  80^pcr2 

7rln  10(1  + 

where  (72  is  the  “volume  parameter”  used  in  the  older  models.  Equation  10  can  be  de¬ 
rived  from  Eq.  IV-69  of  Ref.  [1]  by  taking  73  =  3  and  giving  the  compressional/density 
fluctuation  parameter,  fx,  a  value  of  —1. 
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3.  ACOUSTIC  DEFINITIONS 

This  section  provides  definitions  common  to  the  roughness-  and  volume-scattering 
components  of  the  model. 


3.1  Bistatic  Scattering  Strength 

The  bistatic  scattering  strength  is  the  decibel  equivalent  of  the  bistatic  scattering 
cross  section  per  unit  solid  angle  per  unit  surface  area.  The  scattering  cross  section 
will  be  expressed  as  a  sum  of  terms  corresponding  to  scattering  by  interface  roughness 
and  volume  inhomogeneity;  consequently,  the  scattering  strength  can  be  written  in 
the  form 

s{9s,4>s,0i)  =  I0logio[(^r{9s,(f>s,9i)  +  ay{es,(l)s,9i)]  •  (H) 

The  angular  variables  9s,  (ps,  and  9i  are  defined  in  Fig.  3. 


Figure  3:  Geometry  relevant  to  bistatic  scattering  from  a  transversely  isotropic 
seafloor. 

This  section  uses  the  conventions  defined  in  Ref.  [12]  with  slight  modifications. 
The  grazing  angle  of  the  incident  acoustic  energy  is  denoted  9i,  the  grazing  angle  of 
the  scattered  acoustic  energy  is  denoted  9s,  and  the  change  in  azimuth  due  to  scat¬ 
tering  is  denoted  (ps-  This  angle  will  be  referred  to  as  the  “bistatic  angle.”  In  general, 
two  azimuthal  angles  are  needed,  one  for  the  incident  energy  and  one  for  the  scattered 
energy.  However,  the  model  assumes  that  roughness  is  isotropic  and  volume  inho¬ 
mogeneity  is  transversely  isotropic,  so  that  the  bistatic  scattering  strength  depends 
only  on  the  difference  of  these  two  angles.  Note  that  backscattering  corresponds  to 
the  choice  9s  =  9i,  (ps  =  180°,  and  scattering  in  the  specular  direction  corresponds  to 
9s  —  9i,  (ps  =  0°. 
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3.2  Plane  Waves 


Both  the  roughness  and  volume  scattering  components  of  the  model  employ  first- 
order  solutions  of  the  scattering  problem.  In  several  expressions  to  follow,  complex 
wave  speed  ratios  are  used.  These  are  defined  as 


Up 

l+iip 

(12) 

for  compressional  waves  and 

_ 

1  -f  i6t 

(13) 

for  shear  waves. 

The  solutions  for  plane  waves  propagating  in  the  upward  and  downward  direc¬ 
tions  in  a  homogeneous  medium  can  be  expressed 


exp(ik^  •  r)  , 

where  the  subscript  a  =  p  ov  t  and  denotes  the  wave  vector  for  compressional  or 
shear  waves  in  the  seabed.  The  superscripts  -|-  and  —  denote  the  vertical  component 
of  phase  velocity,  upward  and  downward,  respectively. 

These  three-dimensional  wave  vectors  will  be  decomposed  into  transverse  and 
vertical  components  and  normalized  by  division  by  the  wavenumber  in  water.  That 
is,  if  we  define 

ba  =  ,  (14) 

where 

K  =  —  (15) 

is  the  wavenumber  in  water,  the  normalized  wave  vectors  of  the  incident  downward- 
propagating  compressional  and  shear  waves  in  the  seafloor  are 

b~  =  [cos  6>i,  0,  -PaiOi)]  ,  a  =  p,  t,  (16) 

and  the  normalized  wavevectors  of  the  upward-propagating  scattered  waves  are 

K  =  [cos  Og  cos  4>si  cos  9s  sin  0s,  Pa(^s)]  )  a  =  p,  t.  (17) 

The  vertical  components  of  the  normalized  wave  vectors  in  the  elastic  medium  are 
proportional  to 

Pa{9)  =  cos2(6'),  a  =  p,t  .  (18) 

There  are  two  plane- wave  shear  polarizations,  both  having  particle  displacement 
transverse  to  the  direction  of  propagation.  Only  the  shear  wave  having  particle  dis¬ 
placement  in  a  vertical  plane  is  of  interest,  as  it  is  the  only  one  that  appears  in  the 
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zeroth-order  problem.  The  particle  motion  for  downward-propagating  and  upward- 
propagating  waves  is  in  the  direction  of  the  following  dimensionless  vectors 

b“  =  [Pt(^i),O,cos0i]  ,  (19) 

b+  =  [-Pt{0s)  cos  (j)s,  -Pt{0s)  sin  cos  6>s]  .  (20) 

The  normalization  used  in  Eqs.  (14),  (16),  (17),  (19),  and  (20)  differs  from  that  used 
in  Ref.  [12],  where  the  corresponding  vectors  have  unit  magnitude. 


3.3  Reflection  and  Transmission  Coefficients 


Both  the  the  roughness  and  volume-scattering  portions  of  the  model  require  an 
expression  for  the  reflection  coefficient  for  plane  sound  waves  in  water  incident  upon 
a  flat,  homogeneous  elastic  half-space.  This  reflection  coefficient  is  [19] 


r.(^)  = 


m  - 1 

z{d)  +  1  ’ 


(21) 


where  the  normalized  acoustic  impedance,  z{9),  is 


The  variable  9t  is  the  complex  angle  of  the  shear  wave  in  the  seafloor  material,  de¬ 
pendent  on  the  grazing  angle  in  water,  9: 

cos  29 1  =  2af  cos^  9  —  1  .  (23) 


z{9)  —  ap  sin  9 


cos^  29t  sin^  29t 

nm  m  I 


Equation  22  also  requires 

sin^  29t  =  1  —  cos^  29t  (24) 

as  an  input.  In  later  applications  of  Eqs.  21-24  the  grazing  angle  in  water,  9,  will 
take  on  the  value  of  either  9i  or  9s.  This  will  require  double  subscripting  of  9t,  either 
9ti  or  9ts.  The  reflection  coefficient  is  used  to  define  the  reflection  loss,  L; 

L{9i)  = -20log,M0i)\  .  (25) 


This  is  the  reflection  loss  that  would  be  measured  if  the  seafloor  were  perfectly  flat. 
Seafloor  roughness  will  cause  angular  spreading,  increasing  loss,  and  will  also  cause 
time  spreading. 

The  volume-scattering  portion  of  the  model  also  requires  expressions  for  the 
transmission  coefficients  for  compressional  and  shear  waves.  The  transmission  coeffi¬ 
cient  for  the  scalar  potential  of  compressional  wave  particle  displacement  is  [19] 


rp(0)  = 


COS  29t  sin  9 

PpW 


[1  -  r„(9)l , 


(26) 
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and  the  corresponding  transmission  coefficient  for  shear  is  [19] 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


4.  SMALL-SLOPE  APPROXIMATION 

Yang  and  Broschat  [9]  and  Wurmser  [14]  have  applied  the  small-slope  approx¬ 
imation  to  scattering  by  a  randomly  rough  interface  separating  two  homogeneous 
media,  one  a  fluid  and  the  other  a  lossy  elastic  material.  In  the  present  model,  the 
small-slope  approximation  is  specialized  to  the  isotropic  case.  This  reduces  a  double 
integral  for  the  bistatic  scattering  strength  to  a  single  integral,  and  further  simplifi¬ 
cation  results  from  the  assumption  that  the  coherent  component  of  the  field  is  zero. 
This  is  a  reasonable  approximation  at  high  frequencies  and  was  used  in  the  prede¬ 
cessors  to  this  model  [1].  The  resulting  expression  for  the  bistatic  scattering  cross 
section  is 

OrV.,  4>., «,)  =  4^  l°°  Jo(n)e-'“’"«du  ,  (28) 

where 

P  =  2'^14Q-^  ,  (29) 

where  Jo  is  the  zeroth-order  Bessel  function  of  the  first  kind,  and  Ch  and  a  are 
parameters  of  the  structure  function,  I>(R),  for  two-dimensional  surface  roughness, 
R  =  (x,  y).  For  the  surfaces  characterized  by  the  power-law  spectrum  defined  earlier, 
the  structure  function  is 


D{R)  = 

(30) 

with 

1 

II 

(31) 

and 

2  27r'«;2r(2  -  a) 

^  22“a(l  -  a)r{l  +  a)  ’ 

(32) 

The  variable  Q  is  the  transverse  Bragg  wavenumber,  the  magnitude  of  the  difference 
in  the  transverse  components  of  the  incident  and  scattered  wavevectors. 


Q  =  kw  \J cos^  9s  — 2  cos  Og  cos  di  cos  4)s  +  cos^  9i  -h  IP'  .  (33) 

The  parameter  b  is  assigned  the  value  0.001  to  prevent  numerical  difficulties  that  arise 
if  Q  becomes  too  small.  The  variable  is  the  corresponding  difference  in  the  vertical 
components  of  the  incident  and  scattered  wavevectors 

q3  =  kw{sm6s  +  sin 6*i)  .  (34) 

Equation  28  is  related  to  the  corresponding  Kirchhoff  and  perturbation  expressions. 
The  integral  in  28  is  exactly  the  same  as  that  used  in  the  Kirchhoff  approximation 
in  the  APL-UW  bistatic  model  Eq.  [1],  the  only  difference  being  that  this  integral 
is  to  be  used  for  all  directions  in  the  present  model  while  its  use  was  restricted  to 
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near-specular  directions  in  the  previous  bistatic  model.  It  is  important  to  note  that 
this  integral  is  independent  of  the  properties  of  the  lower  medium;  that  is,  the  same 
integral  is  used  for  the  Dirichlet  (pressme  release  surface),  Neumann  (hard  surface), 
fluid,  and  elastic  problems.  The  coefficient  A  in  Eq.  28  does  depend  on  the  properties 
of  the  lower  medium  and  is  obtained  using  first-order  perturbation  theory  for  the 
corresponding  problem.  Specifically,  the  bistatic  cross  section,  cFrp{9s,<f>siQ%)-,  in  first- 
order  perturbation  theory  is 


<rr,(e„4>.,6i)  =  \A\'‘w(Q) . 


(35) 


The  coefficient.  A,  can  be  obtained  for  the  elastic  case  from  various  sources 
[9, 10, 12, 20, 21].  These  authors  provide  expressions  from  which  A  can  be  obtained 
via  numerical  matrix  inversion.  While  these  expressions  are  readily  coded  using  high- 
level  matrix-oriented  computational  tools,  it  is  possible  to  obtain  the  required  inverse 
analytically,  yielding  the  following  readily  evaluated  scalar  expression  for  A: 

A  =  ^  (Dill  +  r«,(0,)][i  -h  r^(^i)]  +  D2[i  -  r^(^.)][i  +  r^(^i)] 

+£>3[i  +  TMW  -  rt«(0z)]  +  t»4[i  -  r^(^.)][i  -  r^„(^i)])  ,  (36) 

where  rTO(^)  is  the  reflection  coefficient  defined  in  Eq.  21  and  the  variables  Dn  are 

Di  =  —'i-  +  S+— - — - 

a^Op  cos  20 ts  cos  20 a 


of  [(oj  ^  —  2  cos^  Os  — 2  cos^  Oi  -i-  2S)S  2  cos^  Og  cos^  6i 
ttp  cos  29ts  cos  26ti 


D2  = 

£>3  = 


4a^sin^sPt(^5) 
cos  20ts  cos  20ti 


[P^{ei)cosHs  +  {cosHi-S)S] 


? 


4a^  sin  9iPt{9i) 
cos  29ts  cos  29ti 


[P^{9s)cos'^9i  +  {cosHs-S)S] 


5 


(37) 

(38) 

(39) 


£>4  = 


2a^apSm9s  sin 9iPt{9s)Pti9i) 


cos  29ts  cos  29, 


[2{af.  ^  —  2S)S  —  4  cos^  9s  cos^  0i(l  —  2a^a  ^)] 


ti 


—{ttp  —  1)  sin  9s  sin  9i  . 


(40) 


The  variable  S  is 


S  =  cos  9s  cos  9i  cos  0s  .  (41) 

These  expressions  have  not  been  published  elsewhere  but  are  mathematically  equiva¬ 
lent  to  results  given  by  the  authors  cited  above.  The  expression  for  A  is  cast  in  a  form 
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similar  to  that  used  by  Moe  and  Jackson  [15],  in  which  the  flat- interface  reflection 
coefficient,  r„,(^),  appears  prominently.  This  expression  is  obviously  invariant  under 
the  reciprocity  transformation,  in  which  the  incoming  and  outgoing  wave  vectors  are 
interchanged  after  a  sign  reversal.  In  the  present  context,  reciprocity  implies  that  the 
scattering  cross  section  be  unchanged  in  value  by  the  interchange  of  6i  and  Og.  This 
invariance  must  be  obeyed  on  fundamental  grounds  and,  although  it  is  also  obeyed 
by  the  matrix  expressions  given  in  the  references  cited  above,  this  symmetry  is  not 
obvious  upon  inspection  of  any  of  those  references  except  Wurmser  [14]. 
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5.  VOLUME  SCATTERING 


Ivakin  [16]  has  treated  volume  scattering  in  shear-supporting  seafloors  using 
perturbation  theory.  This  approach  has  been  investigated  further  by  Ivakin  and 
Jackson  [11,12]  who  provide  an  expression  for  the  apparent  interface  scattering  cross 
section,  cfvpi^si  due  to  seafloor  volume  scattering  in  the  perturbation  approxi¬ 

mation.  Experience  with  this  expression  has  shown  that  it  sometimes  predicts  unre¬ 
alistically  large  scattering  cross  sections  near  the  specular  direction.  These  values  are 
large  enough  to  violate  energy  conservation  and  signal  a  failure  of  the  perturbation 
approximation.  When  this  occurs,  it  is  reasonable  to  assume  that  volume  scattering 
is  indeed  very  strong,  but  the  actual  scattering  cross  section  will  differ  from  that  given 
by  the  model.  To  accommodate  such  situations,  the  volume  scattering  cross  section, 
^s,0i),  appearing  in  Eq.  11  will  be  limited  to  values  less  than  the  parameter, 
(7o,  as  follows: 

CTviOs,  (t>s,  Oi)  =  [o-"^(0s,  (t>s,  Oi)  +  .  (42) 

In  the  examples  presented  in  this  report,  the  limiting  scattering  cross  section  is  set  to 
the  value  cro  =  0.1,  corresponding  to  a  maximum  scattering  strength  of  —10  dB.  This 
value  was  chosen  to  avoid  violation  of  energy  conservation  in  the  examples  considered 
in  this  report  and  affects  only  the  model  results  for  rough  rock,  rock,  and  cobble  in 
certain  angular  regimes. 

The  perturbation  result  given  by  Jackson  and  Ivakin  [12]  is 


4*st  Oi) 


Im\  Yh  d,,/3d: 


n'p'- 


Vdk'  [(q,  + 
(9.|3  -  Sjs) 


(43) 


where  Im  designates  the  imaginary  part  of  the  complex  quantity  and  is  the 
matrix  of  cross  spectra  for  volume  inhomogeneities  defined  in  Eq.  5.  The  argument 
of  Wppi  can  be  complex,  and  the  resulting  complex  spectrum  values  are  defined  by 
analytic  continuation  of  the  algebraic  expressions  4  and  5. 

The  P  and  /?'  sums  run  over  the  three  types  of  inhomogeneity:  density,  com- 
pressional  wave  speed,  and  shear  wave  speed  (/?,/?'  =  p,p,  t).  The  77  and  rj'  sums  run 
over  the  four  types  of  wave  conversion  caused  by  volume  scattering:  compressional 
to  compressional  (77  =  pp),  shear  to  compressional  {rj  =  pi),  compressional  to  shear 
(77  =  tp),  and  shear  to  shear  (77  =  tt).  Thus,  we  can  put 


77  =  aa', 

where  a  and  a'  run  over  the  two  types  of  waves  (a,  a'  =  p,t).  The  scattering  vectors 
for  the  four  types  of  conversion  are 

Qr,  =  -  b",)  .  (44) 
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These  vectors  give  the  change  in  wave  vector  for  the  corresponding  conversion  and 
appear  as  arguments  in  the  cross  spectra  in  Eq.  43.  The  vertical  components  of 
are  also  used  in  Eq.  43. 

+  Pa'i^i)]  • 


The  coefficients  dr,p  are 
where 


Wr,  =  Ta{0s)ra>i0i)  . 

The  coefficients  Dap  are  elements  of  the  three-column  matrix 

D  =  {Dp\Dp\Dt)  , 


where  Dp  =  D'-\-  Dt/2  and 
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(46) 

(47) 

(48) 


(49) 


(50) 


(51) 


(52) 


where  a,  a'  =  p,  t,  v,  the  normalized  vectors  are  defined  in  Eqs.  16-20,  and  o„  =  0(. 
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6.  MODEL  ILLUSTRATIONS 

In  this  section,  the  model  will  be  compared  with  earlier  backscatter  and  bistatic 
models.  As  the  present  model  has  several  more  parameters,  the  default  parameters 
defined  in  Section  2.4  will  be  used  to  simplify  the  discussion.  The  set  of  parameters 
used  to  illustrate  the  model  is  listed  in  Table  1.  These  examples  can  be  divided  into 
three  classes:  (1)  rough  rock,  rock,  cobble,  and  sandy  gravel,  where  roughness  is  large 
and  shear  effects  are  strong;  (2)  sandy  sediments,  where  roughness  is  moderate  but 
important,  shear  effects  are  small,  and  the  sediments  are  “fast”  (i.e.,  the  sediment 
compressional  wave  speed  exceeds  the  water  sound  speed);  and  (3)  silt,  where  rough¬ 
ness  and  shear  effects  are  negligible  and  the  sediment  is  “slow”  (this  is  also  true  of  clay 
sediments).  In  these  latter  two  classes,  all  parameters  are  assigned  the  default  param¬ 
eters  defined  earher,  except  that  coarse  sand  is  given  nonneghgible  shear  parameters 
taken  from  Ref.  [11].  The  parameters  for  rough  rock  and  rock  are  adapted  from  the 
sedimentary  rock  example  of  Ref.  [11],  for  which  Upp  =  9,  Ot*  =  16,  Opp  =  3,  Utp  =  4, 
and  apt  =  12.  This  choice  gives  perfect  correlation  between  density  and  elastic  mod- 
ulii.  The  aspect  ratio,  a,  is  assigned  the  value  unity,  departing  from  Ref.  [11],  and 
the  spectral  exponents  and  strengths  were  determined  by  trial-and-error  fitting  to  the 
older  backscatter  model.  For  cobble  and  sandy  gravel,  no  guidance  is  available  from 
the  literature,  so  the  defaults  Opp  =  Ou  =  cipp  =  a-tp  =  Upt  =  0,  and  a  =  1  were  used 
along  with  trial-and-error  fitting  to  determine  the  spectral  exponents  and  strengths. 

Because  the  backscatter  model  has  been  built  upon  a  rather  large  data  set  and 
checked  against  data  from  numerous  experiments,  comparison  of  the  predictions  of 
the  new  bistatic  model  with  those  of  the  older  backscatter  model  amounts  to  a  com¬ 
parison  with  data  and  is  a  crucial  test.  This  comparison  is  made  at  a  frequency  of 
20  kHz  for  bottom  types  ranging  from  rough  rock  to  silt.  Comparisons  are  also  given 
between  the  reflection  loss  of  the  fluid  model  [1]  and  the  elastic  model  (Eq.  25) .  Com¬ 
parisons  between  the  older  bistatic  model  [1]  and  the  present  model  are  also  made 
for  incident  grazing  angles  of  10°  and  45°,  at  a  frequency  of  30  kHz.  The  bistatic 
scattering  strength  is  displayed  in  two  ways.  First,  for  a  fixed  incident  grazing  angle, 
the  bistatic  angle  is  set  to  180°,  and  the  scattered  grazing  angle  is  varied  from  0°  to 
180°  (from  backscattering  to  forward  scattering.  Second,  the  scattered  grazing  angle 
is  set  equal  to  a  fixed  value  of  the  incident  grazing  angle,  and  the  bistatic  angle  is 
varied  from  —180°  to  180°  In  the  first  case,  the  backscatter  direction  corresponds  to  a 
scattered  grazing  angle  equal  to  the  incident  grazing  angle  {Og  =  9^,  and  the  specular 
direction  corresponds  to  Og  =  180°  —  6i.  In  the  second  case,  the  backscatter  direction 
corresponds  to  a  bistatic  angle  of  either  180°  or  —180°,  and  the  specular  direction 
corresponds  to  a  bistatic  angle  of  0°. 

All  these  model  comparisons  will  be  discussed  in  order  of  seafloor  type,  beginning 
with  rough  rock  and  ending  with  silt.  The  seafloor  types  were  chosen  purely  for 
illustration  purposes.  In  reality,  each  type  will  consist  of  a  wide  range  of  subtypes 
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whose  parameters  may  differ  significantly  from  those  used  in  this  report.  Thus  the 
results  given  here  should  not  be  considered  representative  of  these  types  of  seafloors 
in  general. 

6.1  Rough  Rock 

Figure  4  compares  the  predictions  of  the  older  backscatter  model  [1]  with  those 
of  the  present  model  at  a  frequency  of  20  kHz.  For  rough  rock,  the  older  backscatter 
prediction  is  purely  empirical,  and  the  approximate  agreement  between  the  older 
and  newer  models  is  a  result  of  the  fitting  process  described  earlier.  As  noted  by 
Tvs, kin  and  Jackson  [11, 12],  a  surprising  result  of  the  inclusion  of  shear  effects  is  that 
volume  scattering  becomes  important  for  rocky  seafloors  and  is  expected,  in  fact,  to 
be  dominant  for  grazing  angles  smaller  than  the  compressional  critical  angle  (64.2°). 
Roughness  scattering  exhibits  a  minimum  for  grazing  angles  somewhat  smaller  than 
the  shear  critical  angle  (40.2°),  while  reflection  loss  (Eq.  25)  reaches  a  peak  in  the 
same  neighborhood.  This  is  a  Rayleigh  wave  effect  [9]  which  provides  a  means  of 
coupling  energy  into  the  rock.  The  volume  scattering  limitation  of  Eq.  (42)  produces 
a  slight  flattening  of  the  angular  dependence  of  scattering  strength. 


Comparison  of  Scattering  Models,  Rough  Rock 


Figure  4:  Comparison  of  predictions  by  the  elastic  and  fluid  scattering  strength  mod¬ 
els  and  reflection  loss  model  for  backscattering  at  20  kHz  from  rough  rock. 

The  old  and  new  reflection  loss  predictions  are  quite  different.  While  the  fluid 
model  predicts  negligible  loss  for  grazing  angles  below  the  compressional  critical  angle. 
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the  elastic  model  predicts  substantial  loss  in  this  region.  This  behavior  is  in  qualitative 
agreement  with  the  conclusions  of  Keenan  et  al.  [22],  although  no  data  are  available 
from  controlled  experiments  to  check  this  prediction. 


Comparison  of  Bistatic  Scattering  Models,  Rough  Rock 


Bistatic  Angle  (deg) 

Figure  5:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering  models 
at  30  kHz  and  a  10°  incident  grazing  angle  for  rough  rock. 


Figures  5  and  6  compare  predictions  of  the  older  and  newer  bistatic  models  for 
incident  grazing  angles  of  10°  and  45°,  respectively.  The  newer  model  shows  a  narrow 
peak  in  the  specular  direction  at  the  smaller  (10°)  incident  grazing  angle  (Fig.  5). 
This  peak  is  due  to  interface  scattering  and  rises  above  a  plateau  owing  to  volume 
scattering.  This  plateau  is  limited  to  a  value  10  log^o  <7o  =  — 10  dB  by  Eq.  42.  While  it 
seems  safe  to  assert  that  volume  scattering  is  important  in  this  region,  the  predicted 
level  and  flat  angular  dependence  seen  in  Fig.  5  axe  suspect  owing  to  the  ad  hoc 
nature  of  Eq.  42. 


6.2  Rock 

The  input  parameters  used  for  rock  are  the  same  as  those  used  for  rough  rock  ex¬ 
cept  that  the  roughness  and  volume  inhomogeneity  levels  are  changed.  The  backscat- 
tering  behavior  (Fig.  7)  and  bistatic  scattering  behavior  (Figs.  8  and  9)  are  similar 
to  those  for  rough  rock,  while  the  reflection  loss  curve  is  identical. 
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Comparison  of  Bistatic  Scattering  Models,  Rough  Rock 
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Figure  6:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering  models 
at  30  kHz  and  a  45°  incident  grazing  angle  for  rough  rock. 


Comparison  of  Scattering  Models,  Rock 


Figure  7:  Comparison  of  predictions  by  the  elastic  and  fluid  scattering  strength  mod¬ 
els  and  reflection  loss  model  for  backscattering  at  20  kHz  from  rock. 
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Comparison  of  Bistatic  Scattering  Models,  Rock 


Figure  8:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering  models 
at  30  kHz  and  a  10°  incident  grazing  angle  for  rock. 


Comparison  of  Bistatic  Scattering  Models,  Rock 
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Figure  9:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering  models 
at  30  kHz  and  a  45°  incident  grazing  angle  for  rock. 
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6.3  Cobble 

The  density  and  compressional  wave  parameters  for  cobble  were  taken  from  Ref. 
[1]  (compare  parameters  in  Figs.  1  and  10),  and  the  shear  and  roughness  parameters 
were  adjusted  to  fit  the  older  backscatter  model  (Fig.  10).  The  resulting  shear  speed 
is  slightly  faster  than  the  water  sound  speed.  Unfortunately,  no  measured  values  are 
available  for  comparison.  Volume  scattering  is  relatively  unimportant  in  this  case, 
and  the  reflection  loss  of  the  new  model  is  much  larger  than  that  of  the  fluid  model 
described  in  Ref.  [1].  The  bistatic  scattering  behavior  (Figs.  11  and  12)  is  similar 
that  for  rock,  with  a  narrow  peak  in  the  specular  direction  rising  above  a  plateau  of 
volume  scattering.  This  plateau  is  narrower  than  in  the  rough  rock  and  rock  examples, 
indicating  that  the  volume  perturbation  approximation  is  providing  believable  results 
over  a  larger  angular  regime. 


Comparison  of  Scattering  Models,  Cobble 
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Comparison  of  Reflection  Loss  Models,  Cobble 


Figure  10:  Comparison  of  predictions  by  the  elastic  and  fluid  scattering  strength 
models  and  reflection  loss  model  for  backscattering  at  20  kHz  from  cobble. 
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Comparison  of  Bistatic  Scattering  Modeis,  Cobble 


Figure  11:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  10°  incident  grazing  angle  for  cobble. 
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Figure  12:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  45°  incident  grazing  angle  for  cobble. 
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6.4  Sandy  Gravel 

For  the  sandy  gravel  example,  all  parameters  except  the  roughness  and  shear 
parameters  were  taken  from  Ref.  [1]  using  the  default  procedure  defined  earlier.  Be¬ 
cause  data  on  shear  waves  in  gravel  were  lacking,  Ut  and  St  were  arbitrarily  assigned 
values  similar  to  those  appropriate  for  coarse  sand  [23-26].  The  two  roughness  pa¬ 
rameters  were  assigned  values  from  the  older  backscatter  model.  Figure  13  shows 
that  the  inclusion  of  shear  increases  the  reflection  loss  somewhat  for  angles  below  the 
compressional  critical  angle  (41.6°).  The  new  model  has  essentially  the  same  bistatic 
behavior  as  the  older  model  of  Ref.  [1]  (Figs.  14  and  15). 


Comparison  of  Scattering  Models,  Sandy  Gravel 
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Figure  13:  Comparison  of  predictions  by  the  elastic  and  fluid  scattering  strength 
models  and  reflection  loss  model  for  backscattering  at  20  kHz  from  sandy  gravel. 
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Comparison  of  Bistatic  Scattering  Models,  Sandy  Gravel 


Figure  14:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  10°  incident  grazing  angle  for  sandy  gravel. 


Comparison  of  Bistatic  Scattering  Models,  Sandy  Gravel 


Figure  15:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  45°  incident  grazing  angle  for  sandy  gravel. 


TM  2-00  27 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


6.5  Coarse  Sand 

The  shear  parameters  for  coarse  sand  were  taken  from  Ivakin  and  Jackson  [11] 
who,  in  turn,  estimated  their  parameters  from  historical  data  [23,  24].  All  other 
parameters  were  taken  by  default  from  Ref.  [1].  As  Figs.  16,  17,  and  18  show,  the 
reflection  and  bistatic  behavior  of  the  fluid  and  elastic  models  is  essentially  the  same. 

The  older  backscatter  model  gives  larger  scattering  strengths  at  small  grazing 
angles,  owing  to  its  use  of  composite  roughness  averaging.  It  has  been  suggested 
that  composite  roughness  averaging  may  be  rather  poor  for  sands  [4],  but  this  issue 
requires  further  examination. 


Comparison  of  Scattering  Models,  Coarse  Sand 


Comparison  of  Reflection  Loss  Models,  Coarse  Sand 
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Figure  16:  Comparison  of  predictions  by  the  elastic  and  fluid  scattering  strength 
models  and  reflection  loss  model  for  backscattering  at  20  kHz  from  coarse  sand. 
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Comparison  of  BistatIc  Scattering  Models,  Coarse  Sand 


Figure  17:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  10°  incident  grazing  angle  for  coarse  sand. 


Comparison  of  Bistatic  Scattering  Models,  Coarse  Sand 


Figure  18:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  Hz  and  a  45°  incident  grazing  angle  for  coarse  sand. 
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6.6  Medium  and  Fine  Sand 

For  medium  and  fine  sand  all  parameters  were  set  by  the  default  procedures, 
which  render  shear  effects  negligible.  Not  surprisingly,  the  predictions  of  the  new 
model  match  those  of  the  older  models  closely  in  all  aspects:  backscatter  and  reflection 
loss  (Figs.  19  and  20)  and  bistatic  scattering  (Figs.  21-24). 

While  volume  scattering  is  relatively  unimportant  for  medium  sand,  it  plays  a 
substantial  role  for  fine  sand.  The  discrepancy  between  backscattering  predictions  by 
the  old  and  new  models  at  small  grazing  angles  is  rather  small  for  fine  sand  owing 
to  the  lower  level  of  roughness  and  the  correspondingly  small  composite  roughness 
averaging  correction. 


Comparison  of  Scattering  Models,  Medium  Sand 
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Figure  19:  Comparison  of  predictions  by  the  elastic  and  fluid  scattering  strength 
models  and  reflection  loss  model  for  backscattering  at  20  kHz  from  medium  sand. 
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Comparison  of  Scattering  Models,  Fine  Sand 


Comparison  of  Reflection  Loss  Models,  Fine  Sand 


Figure  20:  Comparison  of  predictions  by  the  elastic  and  fluid  scattering  strength 
models  and  reflection  loss  model  for  backscattering  at  20  kHz  from  fine  sand. 


Comparison  of  Bistatic  Scattering  Models,  Medium  Sand 


Figure  21:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  10°  incident  grazing  angle  for  medium  sand. 
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Comparison  of  Bistalic  Scattering  Models,  Medium  Sand 


Figure  22:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  45°  incident  grazing  angle  for  medium  sand. 


Comparison  of  Bistatic  Scattering  Models,  Fine  Sand 


Figure  23:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  10°  incident  grazing  angle  for  fine  sand. 


32  TM2-00 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


Comparison  of  Bistatic  Scattering  Models,  Fine  Sand 


Figure  24:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  45°  incident  grazing  angle  for  fine  sand. 

6.7  Silt 

Again,  shear  effects  are  ignored,  and  the  parameters  are  taken  from  Ref.  [1] 
using  the  default  procedures.  This  is  an  example  of  a  “slow”  sediment  for  which 
volume  scattering  dominates  except  near  normal  incidence  for  backscattering  and 
near  the  specular  direction  for  bistatic  scattering.  In  all  cases,  the  elastic  model  is 
nearly  indistinguishable  from  the  older  fluid  models  (Figs.  25-27). 
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Comparison  of  Scattering  Models,  Sill 


Figure  25:  Comparison  of  predictions  by  the  elastic  and  fluid  scattering  strength 
models  and  reflection  loss  model  for  backscattering  at  20  kHz  from  silt. 
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Figure  26:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  10°  incident  grazing  angle  for  silt. 
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Comparison  of  Bistatic  Scattering  Models,  Silt 


Figure  27:  Comparison  of  predictions  by  the  elastic  and  fluid  bistatic  scattering 
models  at  30  kHz  and  a  45°  incident  grazing  angle  for  silt. 
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7.  CONCLUSIONS 

The  elastic  scattering  model  combines  the  best  features  of  three  older  fluid  mod¬ 
els  (for  backscattering,  reflection,  and  bistatic  scattering)  in  a  single  model  and  treats 
a  wider  range  of  seafloor  types.  It  matches  the  older  models  in  the  regime  where  they 
are  well  tested,  but  experimental  tests  are  needed  for  the  extreme  seafloor  types, 
rock  and  very  coarse-grained  material  (cobble  and  gravel).  In  particular,  the  unusual 
prediction  that  volume  scattering  is  important  in  rock  should  be  examined  closely. 
In  addition  to  such  experimental  tests,  mathematical  tests  are  needed  for  these  same 
cases  to  determine  the  accuracy  of  the  small-slope  and  volume  perturbation  approxi¬ 
mations.  It  is  clear  that  the  volume  perturbation  approximation  sometimes  fails  near 
the  specular  direction  for  the  parameters  used  here  for  rough  rock  and  rock. 
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